Duality and the geometric measure of entanglement of general multiqubit W states 
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We find the nearest product states for arbitrary generalised W states of n qubits, and show that the nearest 
product state is essentially unique if the W state is highly entangled. It is specified by a unit vector in Eu- 
clidean n-dimensional space. We use this duality between unit vectors and highly entangled W states to find the 
geometric measure of entanglement of such states. 
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a. Introduction. Quantifying entanglement of multipar- 
tite pure states presents a real challenge to physicists. Inten- 
sive studies are under way and different entanglement mea- 
sures have been proposed over the years [1-6]. However, it is 
generally impossible to calculate their value because the def- 
inition of any multipartite entanglement measure usually in- 
cludes a massive optimization over certain quantum protocols 
or states [7-9]. 

Inextricable difficulties of the optimization are rooted in a 
tangle of different obstacles. First, the number of entangle- 
ment parameters grows exponentially with the number of par- 
ticles involved [10]. Second, in the multipartite setting sev- 
eral inequivalent classes of entanglement exist [11, 12]. Third, 
the geometry of entangled regions of robust states is compli- 
cated [13]. All of these make the usual optimization methods 
ineffective [13-15]. Concise and elegant tools are required to 
overcome this problem. 

A widely used measure for multipartite systems is the ge- 
ometric measure of entanglement E g [16], i.e. the distance 
from the nearest product state. For an n-part pure state tp it is 
defined as E g (tp) = —2\ng(tp), where the maximal product 
overlap g(ip) is given by 



g(tp) = max 

Ml,U 2 ,...,Mfe 



\uiu 2 .-u k ) |, 



and the maximization is performed over all product states. 
The maximal product overlap has many remarkable applica- 
tions. Among them are: it singles out the multipartite states 
applicable for perfect quantum teleportation and superdense 
coding [13], it can create a generalized Schmidt decomposi- 
tion for arbitrary n-part systems [17], it identifies irregular- 
ity in channel capacity additivity [18], it quantifies the diffi- 
culty of distinguishing multipartite quantum states by local 
means [19], it is a good entanglement estimator for quan- 
tum phase transitions in spin models [20], it detects a one- 
parameter family of maximally entangled states [21], and it 
can be easily estimated in experiments [22]. 

In what follows states with g 2 > 1/2 are referred to as 
slightly entangled, states with g 2 < 1/2 are referred to as 



highly entangled and states with g 2 = 1/2 are referred to as 
shared quantum states. In this Letter we show how to calculate 
the maximal product overlap of an arbitrary W state [11]. The 
method is to establish a one-to-one correspondence between 
highly entangled W states and their nearest product states. 
Consider first generalized Greenberger-Horne-Zeilinger 

states [23], i.e. states that can be written |GHZ) = a|0 0) + 

6| 1 ... 1) in some product basis. Such states are fragile under 
local decoherence, i.e. they become disentangled by the loss 
of any one party, and they are not highly entangled in the sense 
defined above. The geometric measure of these states is com- 
puted easily since the maximal overlap simply takes the value 
of the modulus of the larger coefficient, \a\ or \b\ [24]. Ac- 
cordingly, the nearest separable state is the product state with 
the larger coefficient. Thus many generalized GHZ states with 
different maximal overlaps can have the same nearest product 
state. 

Consider now generalized W-states [25], which can be writ- 
ten 

|W„) = ci|100...0) + c 2 |010...0) + • • • + c„|00...01). (1) 

Without loss of generality we consider only the case of pos- 
itive parameters c k since the phases of the coefficients Cfe 
can be eliminated by redefinitions of local states \lk), k = 
1, 2, n. The states (1) are robust against decoherence [26], 
i.e. loss of any n — 2 parties still leaves them in a bipar- 
tite entangled state. Surprisingly, if the state is slightly en- 
tangled, then we have the same situation as for generalized 
GHZ states: the maximal overlap is the largest coefficient 
and, as before, many states can have the same nearest prod- 
uct state [27]. However, the situation is changed drastically 
when the state is highly entangled. The calculation of the 
maximal overlap in this case is a very difficult problem and 
the maximization has been performed only for relatively sim- 
ple systems [9, 14, 16, 24, 27-30]. 

On the other hand, different highly entangled W-states have 
different nearest product states. This makes it possible to 
map the W-state to its nearest product state and quickly ob- 
tain its geometric measure of entanglement. More precisely, 
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we construct two bijections. The first one creates a map be- 
tween highly entangled n-qubit W states and n-dimensional 
unit vectors x. The second one does the same between n- 
dimensional unit vectors and n-part product states. Thus we 
obtain a double map, or duality, as follows 

|W„) <-» x o \u{) <8 \u 2 ) ® • • ■ ® \u n ). (2) 

The main advantage of the map is that if one knows any of 
the three vectors, then one instantly finds the other two. 

b. Classifying map. Now we prove a theorem that pro- 
vides a basis for the map. 

Theorem 1. Let |W„) be an arbitrary W state (1) with non- 
negative real coefficients c it and let \u\) ® \u 2 ) ® ■ ■ ■ <B> \u n ) 
be its nearest product state. Then the phase of \uk) can be 
chosen so that 

\u k ) = sm0 k \O) + cos0 fe |l), < k < |, k = 1, 2, ...,n. 
where 

cos 2 6»i + cos 2 2 H h cos 2 9 n = 1. (3) 

Proof. The nearest product state is a stationary point for the 
overlap with |W n ), so the states \u k ) satisfy the nonlinear 
eigenvalue equations [9, 16, 17] 

{u\u 2 ■■ - Uk- • -u„|W„) = g\u k ); k = 1,2, • • • ,n (4) 

where the caret means exclusion. We can choose the phase of 
\u k ) so that \u k ) = sin6> fe |0) + e^* cos6> fc |l), and then (4) 
gives the pair of equations 

c k Y[ sin 6j = .9 e ^ fc cos &k, (5a) 



e~ l<t " ci cos 0i Y[ sin 0j = 9 sin k . (5b) 

Eq. (5a) shows that ge l ^ k is real, so ^ = — arg(g) 
is independent of k. Then the modulus of the overlap 
| (ui ■ ■ ■ u„|W„) | is independent of <j>, so we can assume that 
4> = 0. Now multiplying eq.(5b) by sm0 k and using eq.(5a) 
gives Eq.(3). □ 

Thus the angles cos^ define a unit n-dimensional Eu- 
clidean vector x. We can also define a length r as follows. 
From Eq.(5a) it follows that the ratio s'm20 k /c k does not de- 
pend on k. If this ratio is non-zero we can define 

1 _ sin26>i sin26> 2 sin26>„ 

- = = = ■ ■ • = . (6) 

r ci c 2 c„ 



c. Highly entangled W states. Equations (5) admit a 
trivial solution sin20 k = 0, k = 1, 2, • • • , n and a special 
solution with nonzero values of all sines. The trivial solution 
gives the largest coefficient of |W„) for the maximal overlap 
and is valid for slightly entangled states. We consider them 
later and now focus on the special solutions. From Eq.(6) it 
follows that 



cos 2 9 k = - |l±yi--|J , fc = 1,2,.-- ,n. (7) 

The plus sign means that cos2#fc > 0. Then from Eq.(3) it 
follows that this is possible for at most one angle; specifically, 
we prove that if cos 20 k > for some k, then c k is the largest 
coefficient in Eq.(l). Suppose cos26>fc > but c k is not the 
largest coefficient and there exists a greater coefficient, say c;. 
Then from Eq.(6) it follows that sin 26*/ > sin 20 k > and 
consequently | cos 20 1 \ < |cos20fc|. Now we rewrite Eq. (3) 
as follows: 

- cos 26<i - cos 26*2 cos 20 n = n - 2. (8) 

From | cos 20i \ < | cos 20 k | and cos 20 k > it follows that 
— cos 20 k — cos 20i < which is in contradiction with Eq.(8). 
Thus Cfe must be the largest coefficient. 

Without loss of generality we assume that < c\ < • • • < 
c„. Then in (7) we must take the — sign for k = 1, . . . , n — 1 
and (3) becomes 

We will denote the left-hand sides of these equations as f±(r). 
We also use fo(r) to denote this expression without the last 
term. The function r(ci, C2, c n ) defined by /+(?*) = n — 2 
is a completely symmetric function of the state parameters c k . 
In contrast, the function defined by /_ (r) = n — 2 is an asym- 
metric function since its dependence on the maximal coeffi- 
cient c n is different. Thus in equation (9) the upper and lower 
signs describe symmetric and asymmetric entangled regions 
of highly entangled states, respectively. 

For highly entangled states, eqs. (9)± uniquely define r as 
a function of the state parameters c k . More precisely, 

Theorem 2. There are two critical values r\ and r 2 of the 
largest coefficient c n , i.e. functions o/ci, . . . , c„_i such that 

I- If Cn < Ti> there is a unique solution of (9 + ) and no 
solution of (9-); 

2. If c n = r\, both (9+) and (9~) have a unique solution, 
the same for both; 

3. If n < c n < r% there is no solution of (9)+ and a 
unique solution of (9-); 

4. If c n > r 2 , neither (9+) nor (9-) has a solution. In this 
case the state | W„) is slightly entangled. 
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The value n is the solution of /o(?*i) = ti — 2, which exists 
and is unique since / (c n _i) < n — 2 and /o(r) — > n — 1 
monotonically as r — >• oo; and r2 is defined by 



= cf + 



(10) 



Then ^2 > r\, for f 0(1*2) > ?i — 2 = /o(?"i) using -y/a; > x for 
< x < 1. Since /o is an increasing function of r, it follows 
that r2 > ri . Now the theorem follows from the following 
properties of the functions f±(r)(f'_ is the derivative of /_): 

1 . fo and /+ are monotonically increasing functions of r. 

2. ,f+(r) — > n as r — > 00. 

3. If c„ < n, /+(c„) = / (c„) < /o(n) = n - 2. 

4. If c n > n, then /+(r) > n — 2 for all r > 7*1. 

5. If c n < n, then /_(c„) < n — 2. 

6. If c n > n, then /_(c„) > n — 2. 

7. If c n < T2, then /_ (r) < n — 2 for large r. 

8. If c„ > r2 then /_ (r) > n — 2 for large r. 

9. /i(c„ + e) <0 for small e. 

10. If c„ > r 2 , then /i(r) < for all r > c n . 
These properties are illustrated in Figure 1. 
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FIG. 1 : (Color online) The behaviour of the functions f± for five- 
qubit W states. The function /+ (r) (dotted line) and /_ (r) (solid 
line) are plotted against r in the four cases c„ < ri, c n = ri, n < 

Cn < r 2 and c„ = r 2 . 



t/. Geometric measure. We can now identify the nearest 
product state, and the largest product state overlap g(\W n }), 
for any W-state |W„), as follows. 

Theorem 3. If c n > 1/2, the state |W„) defined by (1) is 
slightly entangled. Its nearest product state is |0 . . . 01), with 
overlap ff(|W n )) = c„. 

c„ < 1/2, the state |W„) is highly entangled and has 
nearest product state 



\ Ul )...\u n ) where \u k ) = sin0 fc |O) + |e l * cos6» fc |l) 
w/f/z which its overlap is 

g = 2r sin 0i sin #2 • ■ • sin 9 n . 



(11) 
(12) 



i/ere r is the solution of(9)±, whose existence and uniqueness 
are guaranteed by Theorem 2; the phase <f> is arbitrary; and 
6k is given by (7) with the — sign for k = 1, . . . , n — 1, the — 
sign for k = n ifr satisfies (9+), the + sign if r satisfies (9_). 

Proof. The nonlinear eigenvalue equations (4) always have n 
solutions 



.9 = Cfc, 



|0) if f ^ k, 
II) if i = k 



If c„ > 2, i.e. in case (4) of Theorem 2, there are no other 
stationary values, so the largest overlap <?(|W„)) equals the 
largest coefficient c„, the corresponding product state being 
I0...01). 

If Cn < 1/2 there is another stationary value given by (12). 
We will now show that this is larger than any of the triv- 
ial stationary values Cfc. We use the following inequality: If 
yi, . . . ,y n are real numbers lying between and 1, and satis- 
fying 7/1 H h y n < 1, then 



(l-2/l)(l-y2)" ■ (l-2/n) > 1-V1-V2- 



Vn- (13) 



This is readily proved by induction. We can apply ( 1 3) to n — 1 
terms of Eq.(3) to get 



(1-cos 2 6>i) • • • (1-cos 2 6 n -i) > 1-cos 2 0i- 



-cos 



or 



sin 0i sin 02 sin 0„_i > cos 



(14) 



Now from Eq.(5a) it follows that g 2 > c n . Thus g is the 
maximal product overlap, and the nearest product state is 
\ui) ■ ..u n ). 

Next we prove that if |W„) is normalised, then g 2 < 1/2. 
For this we need another inequality: If y\ , . . . , y n are real 
numbers lying between and 1, and satisfying yx + ■ ■ ■ +y n = 
n — 1, then 



2/1 + — h y n > yl h \-yl + 2 ym 



This can also be proved by induction. 



From (6), and using c\ 



3? = 1, we find 



sin 2 20i 



sin 2 20„ 



Hence (12) gives 



2/12/2 ■ ■ ■ 2/n 



2/l(l — 2/l) H 2/„(l - yn) 



(15) 



(16) 



(17) 



where j/^ = sin 0^. But yi 



2/n = J7—1, so the 



inequality (15) applies, and gives g < 1/2. 



□ 



Finally, we summarise the correspondence between highly 
entangled W-states, their nearest product states, and unit vec- 
tors in t™. 

Theorem 4. There is a 1:1 correspondence between highly 
entangled states |W„) defined by (1), their nearest product 
states with real non-negative coefficients, and unit vectors x £ 
R n with 0<x k < 1/V2 (k = 1, . . . ,n - 1), < x n < 1. 
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Proof. By Theorem 3, |W n ) is highly entangled if and only 
if c n < 1/2. If this is the case, Theorem 1 and (7) show 
that its nearest product state is of the form (11) where x = 
(cos 9i , . . . , cos n ) is a unit vector in W 1 in the region stated. 
The angles Ok are given in terms of the coefficients Cft by 
(6), in which r is a function of the coefficients which, by 
Theorem 2, is uniquely defined. The nearest product states 
|ui)|u2) • ■ • \u n ) are determined by these angles, up to a phase 
<f>, by \uk) = sin(9fc|0) + c 4 * cos0fc|l), so there is only one 
nearest product state with real non-negative coefficients, and 
only one unit vector x, for each highly entangled state |W„). 
Conversely, given a unit vector x = (cos 6\, . . . , cos 6 n ), the 
quantity r is determined by (16), and then the coefficients 
ci , . . . , c n are determined by (6). Thus the correspondences 
(2) are bijections. □ 

The equations (9±) cannot always be explicitly solved to 
give analytic expressions for r in terms of the coefficients Cfc. 
However, in some cases, including all states for n = 3, ex- 
plicit solutions can be obtained. Then the angles Of. can be 
calculated from (6) and eq.(12) gives a formula for the max- 
imal product overlap g(\W n ))- This formula is valid unless 
any of the angles 8^ vanishes, and restores all known results 
for the maximal overlap of highly entangled W states. When 
n = 3 it coincides with the formula (31) in Ref.[9]. When 
Ci = C2 = • • • = c n it coincides with the formula (52) in 
Ref . [24] . And when n = 4 and C3 = C4 it coincides with the 
formula (37) derived in Ref. [27]. 

When max(cj, c|, ■ • • c£) = r| = 1/2 the two expressions 
for <?(|W n }) given in Theorem 3 coincide; these states are 
shared quantum states. The nearest product states and max- 
imal overlaps of shared states are given by the first case of 
Theorem 3, but also they appear as asymptotic limits of the 
second case. Indeed, at the limit 9 n — > we have 

lim 2rsin#„ — > c„, lim 2rcos#/ c — ¥ cu, k 7^ n. (18) 



Thus the angle 6 n vanishes and the length of the vector r 
goes to infinity, but their product has a finite limit. Substi- 
tuting these limits into Eq.(3) one obtains c n — > r|. Therefore 
entangled regions of highly and slightly entangled states are 
separated by the surface c n = 1/2; for states on the surface, 
r —> 00. All of these states can be used as a quantum channel 
for the perfect teleportation and superdense coding [13]. 

e. Summary. We have constructed correspondences be- 
tween W states, n-dimensional unit vectors and separable 
pure states. The map reveals two critical values for quantum 
state parameters. The first critical value separates symmetric 
and asymmetric entangled regions of highly entangled states, 
whiles the second one separates highly and slightly entangled 
states. The method gives an explicit expressions for the ge- 
ometric measure when the state allows analytical solutions, 
otherwise it expresses the entanglement as an implicit func- 
tion of state parameters. 

It should be noted that the bijection between W states and 
?i-dimensional unit vectors is not related directly to the geo- 
metric measure of entanglement. Therefore it is possible to 
extend the method to other entanglement measures. To this 
end one creates an appropriate bijection between unit vec- 
tors and optimization points of an entanglement measure one 
wants to compute. This work is in progress. 
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